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Abstract 

We study the Polyakov loop correlator in the (2+1) dimensional Z2 gauge model. 
An algorithm that we have presented recently, allows us to reach high precision results 
for a large range of distances and temperatures, giving us the opportunity to test 
predictions of the effective Nambu-Goto string model. Here we focus on the regime 
of low temperatures and small distances. In contrast to the high temperature, large 
distance regime, we find that our numerical results are not well described by the two 
loop-prediction of the Nambu-Goto model. In addition we compare our data with 
those for the SU(2) and SU(3) gauge models in (2+1) dimensions obtained by other 
authors. We generalize the result of Liischer and Weisz for a boundary term in the 
interquark potential to the finite temperature case. 



Recently a renewed interest has been attracted by the effective string description of the 
interquark potential in lattice gauge theories (LGT) [I]-[TT]. On one hand, in E] high 
precision simulations of the interquark potential were run, using algorithms jU IH Ej that 
allow to measure the Polyakov loop correlation function with highly reduced variance. On 
the other hand, in [H] the lowest excited states in the spectrum of the confining flux tube 
in the (2+1) dimensional Z 2 gauge model were studied as well, using a variant of earlier 
methods. These works enabled one to test some longstanding conjectures on the effective 
string description of the interquark potential, and also triggered some new theoretical effort 
toward a better understanding of the effective string model itself. While it seems by now 
understood that the first order correction to the interquark potential is given by the so 
called Liischer term ^2] > several other issues are still open and require further investigation. 
A tentative list of these open problems includes: 

• The determination of higher order terms (and possibly the full functional 
form) of the effective string action. 

The presence of the Liischer term in the interquark potential only tells us that, at 
leading order, the effective string action is simply a two-dimensional quantum field 
theory of d — 2 free bosonic fields (one for each transverse degree of freedom of the 
fluctuating string). It does not help to identify the higher order terms in the action, 
which should describe the string self-interaction. To this end, one has to evaluate 
higher order corrections (i.e. higher order powers in 1/R, R being the interquark 
distance) in the interquark potential. We addressed this problem in two recent papers 
[3JE], testing the simplest possible effective action, i.e. the Nambu-Goto string (see 
below for a detailed discussion) in the large distance, finite temperature regime of the 
3d Ising gauge model. We found an impressive agreement between Monte Carlo data 
and the prediction from the Nambu-Goto action truncated at second order l , but 
apparently no room was left for higher order corrections, which should necessarily be 
present if the Nambu-Goto string is the correct picture. 

• The universality issue. 

As we mentioned above, the first order correction to the interquark potential (i.e. 
the Liischer term) shows an impressive degree of universality and has been detected 
with the precise numerical value predicted by the effective string theory in all of the 
models studied up to now (ranging from the 3d gauge Ising model to the SU(2) and 
SU(3) models, both in d=3 and in d=4). Which is the situation for higher order 
corrections? Do they also show the same universal behaviour? Preliminary results 
suggest that this is not the case |2|, however these evidences were obtained in the 
short distance regime only (at present only in the case of the 3d Ising model one can 
reach the large distance regime), where non- universal boundary terms are present 
(see last item below) and make the analysis much more difficult. 

1 A similar agreement was observed in where the large distance regime of the interface potential in 
the 3d Ising model was studied 
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• The boundary corrections. 

As it was recently pointed out by Liischer and Weisz, the presence of two boundaries 
(the two Polyakov loops) in the finite temperature geometry makes it necessary to 
include in the possible higher order terms of the string action also contributions lo- 
calized at these boundary (which we shall call in the following "boundary terms") 
which are proportional to a non-universal parameter, which we shall denote as b in 
the following. These terms induce corrections to the expected form of the interquark 
potential (which we shall call "boundary corrections" in the following ): these cor- 
rections are proportional to b. The open problem in this case is to evaluate b in the 
various LGT's, to check if it can be compatible with zero and more generally if it 
depends on some relevant feature of the underlying LGT. To this end, it would be 
very useful to derive the whole functional form (i.e. the explicit dependence from R 
and L) of the boundary correction. 

• High precision test of the Liischer correction. 

While, as we mentioned above, it is by now accepted that the first order correction 
to the interquark potential is given by the Liischer term, a high precision test of this 
result in the finite temperature regime is still missing. It would be interesting to 
have a quantitative estimate both of the range of values of the interquark distance in 
which the Liischer term well describes the interquark potential and of the numerical 
uncertainty in the determination of its value. This is particularly important in view 
of the recent results j2 E! which show that the bosonic string seems to work at 
surprisingly short distances. 

• A "Casimir energy paradox"? 

Finally, we mention the fact that in jH] the authors presented numerical results show- 
ing that in a distance range (below 1 fm) where the Liischer term in the interquark 
ground state potential is already observed, they only found a few stable excitations 
of the confining flux, whose spectrum turned out to be grossly distorted compared 
with the string prediction. 

The aim of this paper is to address some of these issues by looking once again at the 3d 
gauge Ising model, which allows to perform very precise simulations while keeping all of the 
possible sources of systematic error under control. To this end, we performed a set of new 
simulations in the short distance, low temperature regime of the interquark potential, 
which we used first to make a high precision test of the Liischer term (as discussed above) 
and then to study the higher order terms of the effective string action. In this respect 
the present paper can be considered as the completion of our previous paper [3] in which 
we performed a similar study in the large distance, high temperature regime of the 
interquark potential. Our results can also be compared with the findings of an analogous 
study, recently published in 0. A preliminary account of our results can be found in [TI] . 
A non-trivial problem one has to face when looking at the higher order effective string 
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terms is the presence of the non-universal boundary correction discussed above. In order 
to deal with this problem we evaluated the functional form of the boundary correction (a 
result which is rather interesting in itself as we shall see below), in the framework of the 
zeta-function regularization. In this way, we may predict the large distance behaviour of 
the boundary correction, and use the high precision results of |4. to fix the non-universal 
constant in front of it. As we shall see, our data strongly suggest that this value is 
compatible with zero for the 3d gauge Ising model. 

A side consequence of this analysis is that we shall be able to show in a rather unam- 
biguous way that, at least in the short distance regime the interquark potential, in the 3d 
Ising model cannot be described by a simple Nambu-Goto string. We shall then conclude 
the paper with a tentative comparison of our results with the short distance behaviour of 
other LGT's trying to shed some light on the universality issue discussed above. 

This paper is organized as follows. We start in sect, ^with a short summary of the 
effective string description of the interquark potential. Then in the following two sections 
we compare our new Monte Carlo results with the effective string predictions in the short 
distance regime. In sect. El we discuss the first order correction (the "Liischer" term), 
while in sect. El we address the issue of the higher order corrections. Sect. El is then devoted 
to the study of the functional form of the boundary correction. These results are used in 
sect. El to study the boundary correction in the large distance regime. Sect. EH is devoted 
to some concluding remark and to the comparison of our Ising results with those obtained 
in the SU(2) and SU(3) gauge theories in d = 3. Finally, the detailed derivation of the 
functional form used in sect. E]is reported in the Appendix. 

1 Summary of known theoretical results on the effec- 
tive string description of the interquark potential 

We refer to our previous papers [3 El for a detailed discussion about the effective string 
picture and its peculiar realization in the finite temperature geometry (i.e. in the case 
of Polyakov loop correlators). Here, we only recall some basic facts (more details can be 
found in the appendix). 

• In finite temperature LGT's the interquark potential can be extracted by looking at 
correlators of Polyakov loops in the confined phase. The correlator of two loops P(x) 
at a distance R and at a temperature T = 1/L [L being the size of the lattice in the 
compactified "time" direction) is given by 



The free energy F(R, L) is expected to be described, as a first approximation, by the 
so called "area law" 



G(R) = (P{x)P\x + R)) = e" F(i? ' L) 



(1) 



F(R,L) ~ crLR + k(L) 



(2) 
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where a denotes the string tension and k(L) is a non- universal constant depending 
only on L. 



According to the effective string picture, in the rough phase of the theory one must 
add to the area law of eq. (J2J) a correction due to quantum fluctuations of the flux tube 
("effective string corrections"). Such a correction is expected to be a complicated 
function of R and L, but if one neglects the contributions due to the string self- 
interaction terms and from the boundary terms which we shall discuss below, then 
one finds 

F(R, L) ~ F q (R, L) = aLR + k(L) + F* (R, L) (3) 

with 

F 9 1 ( J R,L) = (rf-2)log(r / (r)) ; -ir = ±= , (4) 



where (d — 2) is the number of transverse dimensions, rj denotes the Dedekind eta 
function 



tj(t) = q 



q" 



q 



(5) 



n=l 



and the labels q and 1 in F^ recall that this is the first order term in the expansion 
of the quantum fluctuations of the flux tube. Eq. (JHJ) is referred to as the "free 
string approximation" . 

The Dedekind function has two different expansions — both of which turn out to be 
convergent for any finite value of L and R — which are most suitable for the two 
regimes 2R < L and 2R > L, respectively. These expansions are related to each 
other by the so called "modular transformation" . 



For 2R< L 



F}(R,L) 



ttL 
2AR 



En 



71=1 



-TrnL/R\ 



(d-2) 



The first term of this expansion is the well known "Luscher term" 
For 2R> L 



F}{R,L) 



7TR 

QL 



71=1 



(d 



(6) 



(7) 



In this case, the first term of the expansion is proportional to R and acts to 
lower the string tension. 

Unless we are in the intermediate region R ~ L/2, the exponentially decreasing 
terms which appear in eq.s © and (JJJ) can be neglected. We shall largely use this 
approximation in the following. 
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At large enough temperatures, i.e. for small values of L, or small enough values 
of R, higher order terms in the expansion of the flux tube quantum fluctuations 
become important and cannot be neglected. These terms encode the string self- 
interaction and depend on the particular choice of the effective string action. The 
simplest proposal (discussed in [SI El) is the Nambu-Goto action, in which the string 
configurations are simply weighted by the area of the world-sheet. Its contribution 
to the free energy turns out to be (fixing for simplicity d — 3) (see ^3] and [3]) 

where E 2 and are the Eisenstein functions. The latter can be expressed in power 
series 

oo 

E 2 {t) = l-24^aH<f (9) 

n=l 

oo 

Ei(r) = l + 240^<7 3 (n)g" (10) 

n=l 

q = e 2 ™ , (11) 

where cr(n) and 03 (n) are, respectively, the sum of all divisors of n (including 1 and 
n), and the sum of their cubes. 



2 High precision test of the Liischer term 
2.1 The simulations 

In order to perform a high precision test of the Liischer term we run a set of Monte Carlo 
simulations in the short distance regime. In tab. Q we report some information on our 
simulations. The choice of very large lattice sizes in the "temperature" direction (i.e. 
the direction of the Polyakov loops) ensures us that we are in the very low temperature 
domain. The temperature of our present simulations ranges from T/T c = 1/20 for L = 80 
at /3 = 0.73107 up to T/T c = 1/5 for L = 20 at (3 = 0.73107. The values for the l/L 
corrections estimated in jlj indicate that for T < the contribution of possible higher 
order effective string terms due to the finiteness of the lattice size in the temperature 
direction is completely negligible within the precision of our data while the data at T = T c /5 
are at the border of our precision. Thus, at least for T < we may neglect possible l/L 
corrections and we can concentrate on the 1/R corrections only. Notice that this is exactly 
the opposite situation with respect to our previous paper jl] in which the temperature was 
much higher and the l/L corrections were dominating. The sample at the smallest value 
of j3: j3 = 0.65608 is characterized by a rather small correlation length and we should 
therefore expect large scaling violations. 
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For a detailed discussion of the algorithm that we used, of the simulation setting and 
more generally on the 3d gauge Ising model, we refer the reader to our previous paper jl]. 
The results of our simulations are reported in tab. Eland in tab. El 2 . 



p 


L c = 1/T C 


L 


N s 




a 


R c 


0.65608 


2 


20 


32 


0.644(1) 


0.20487(1) 


2.71 


0.73107 


4 


20,40,80 


64 


1.456(3) 


0.0440(3) 


5.84 


0.75180 


8 


80 


128 


3.09(1) 


0.010560(18) 


11.92 



Table 1: A few information on our simulations. In the first column the value of (3 is given 
and in the second the inverse of the critical temperature. The third and fourth columns 
contain the values of lattice size in the "temperature" direction L and in the two spacelike 
directions N s . In the last three columns the values of the exponential correlation length 
(obtained by interpolating the values reported in \17\j and \18lj and by the analysis of the 
low temperature series llty). the value of the string tension in dimensionless units, and the 
value of R c defined as R c = a/I-5/ct (see sect. 2.7 of ref. are provided. 



2.2 Test of the Liischer term 

We performed our analysis in three steps. First we made, following what is usually done in 
the literature, a rather naive test by directly fitting the data with a 1/R correction. As we 
shall see below, this choice is correct only in a very narrow range of interquark distances, 
since for large enough values of R the exponentially decreasing terms exp(— irnL/R) con- 
tained in eq. (JBJ) start to matter and cannot be neglected. In the second level of analysis we 
fit the data with the whole functional form of the free bosonic correction, eq. ©, finding 
a remarkable agreement between the data and our theoretical expectations. Finally in the 
third level of analysis we combine the data in such a way that the string tension is elimi- 
nated as parameter. Thus, we can make an absolute comparison, with no free parameter 
to fit, between our data and the effective string prediction. Before going into the details 
of the analysis let us stress a point which will be important in the following. Due to the 
algorithm we used, each number in tab.s 121 and |3] corresponds to a different simulation. 
Hence the sets of data we fit are completely uncorrelated and we can safely trust both on 
the x 2 obtained from the fit and on the best fit results for the parameters. 

2 Notice that also the result at R = 80 listed in tab. |3] and those at R = 40 listed in tab. |3 are 
meaningful (even if 80 > 128/2 and 40 > 64/2) due to the non-trivial mapping of the boundary conditions 
from the Ising gauge to the Ising spin model under duality transformation (see |16| . Sect. 4.3 for a detailed 
discussion of this point). In particular, performing our simulations in the 3d Ising spin model (as we did), 
we neglected the anti-periodic boundary conditions, which would produce the "echo" contribution due to 
the distances (128-80) and (64-40) respectively. 
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Table 2: Values of the ratio between two successive Polyakov loop correlators G(R + 
1)/G(R) for various values of R and L, at (3 = 0.73107 . 

2.2.1 Fitting with a 1/R correction 

We fitted the data reported in tab. |2]and tab. El according to the law 

With this normalization the bosonic string model predicts a = a and b = n/24 = 0.13089... 

If we fit, for a given choice of (3 and L, all the values of R listed in the tables, the 
reduced \ 2 turns out to be very high and the best fit values of a and b are very far 
from the expected values. This indicates that at short R we have higher order (string) 
corrections which are proportional to higher powers of 1/R and compromise the fit. The 
standard way to deal with this type of behaviour is to repeat the fit eliminating the data 
one after the other starting from those with the lowest values of R, until for some value of 
R m in an acceptable (i.e. order unity) reduced \ 2 is finally reached. However it turns out 
that for all our samples, except the one at T = T c /20 this scenario never occurs. Instead 
we see that after a rapid decrease, the Xr starts again to increase as R m in is increased and 
the minimum of this shape never goes below 1. As an example, we report in tab. El the 
results of the fit for the data at f3 = 0.73107 and L = 40. The reason for this behaviour 
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Table 3: Same as tab|2 but for the data at (3 = 0.75180 . 

is that corrections for large values of R are present in the data. These corrections are 
nothing but the exponentially decreasing terms of eq. ©• If we neglect them in the fit, 
we shall never be able to extract the correct value of the Liischer term, as the best fit 
values reported in tab. El clearly show. If the temperature is low enough, there is a range 
of values of R in which both the small R and the large R corrections become negligible 
and the expected value of the Liischer term can be recovered. This range can be found by 
iteratively eliminating the smallest and the largest values of R in the fit and looking for 
a stable and acceptable value of Xr- 

However, if T is not small enough, this procedure does not converge. This is the case 
for instance of our data at /3 = 0.73107 and L = 20. We report in tab.|Hlthe best fit results 
obtained in this way for the other three samples. 

We insisted in this type of fit because it is widely used in the literature and, as we have 
seen, it can easily lead to misleading conclusions. In particular one should be aware that: 

• if L is not large enough, i.e. if the temperature is not small enough, the range of 
interquark distances within which this type of fit can be performed vanishes. The 
precise threshold in temperature depends on the precision of the data. In our case 
the threshold was around T ~ T c /10. 

• If one nevertheless performs the fit outside the allowed region the reduced x 2 becomes 
larger than 1 and the values of the Liischer term which one extracts from the fit turns 
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20 


0.016643(6) 



Table 4: Same as tab|2 but for the data at [3 = 0.65608 . 

out to be (erroneously) different from the expected 7r/24 factor. 

• If one enforces the expected 7r/24 coefficient and performs the fit outside the allowed 
range then the values of the string tension extracted from the fit may show relevant 
systematic deviations from the correct value. 

We think that these observations should be carefully taken into account when dealing 
with the interquark potential extracted from Polyakov loops in any LGT, i.e. also in 
the physically more interesting SU(2) and SU(3) cases, if the precision reached by the 
simulations is high enough. 

2.2.2 Fitting with the whole functional form of the bosonic string correction 

In order to keep into account the next to leading terms in eq. (0) we fitted the data reported 
in tab. Eland tab. El according to the law 

=a+^(F^R+l,L)-F^R,L)) . (13) 

With this normalization a is again the string tension while c = 1 means that we have a 
perfect agreement between the whole bosonic string prediction (Luscher term plus expo- 
nentially decreasing corrections) and the data. In tab.0we report as an example the result 
of the fits in the particular case j3 = 0.75180, L = 80. Similar results are obtained with 
the other values of (3 and L. 
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Table 5: Results of the fit according to eq. ifT^) for the data at (3 = 0.73107 and L = 40. In 
the first column the minimum value of R included in the fit is given. The second column 
contains the reduced v 2 . In the last two columns the best fit results for a and for the 
Luscher term are given, respectively. 

We find a very good agreement between our data and the bosonic effective string 
prediction, keeping in the fit all the interquark distances up to R = 80. The uncertainty in 
our determination of c is less than 1%. The decrease in the reduced x 2 indicates that the 
higher order string contributions can be neglected for R > 20. Looking at the string tension 
results, we see that the best fit value changes by a rather large amount (with respect to 
the quoted errors) as R m i n increases. This observation has two relevant implications: 

• As we mentioned above, a large systematic deviation in the best fit estimates of the 
string tension has to be expected, if the effective string corrections are not properly 
taken into account. 

• Even if the first order effective string correction is properly taken into account, sys- 
tematic deviations in the string tension may still appear if too small values of R are 
included in the fits. 

Our results for the other values of (3 are summarized in tab. |H1 Looking at this table one 
can see that the data do not show scaling violations and that the expected value c = 1 
is found even for (3 = 0.65608 (for which instead any bulk observable shows very large 
scaling violations) ! This represents a further nontrivial test of the whole picture, since the 
effective string picture is indeed expected to hold in the whole range (3 r < [3 < (3 C {(3 T and 
(3 C being respectively the roughening and the deconfinement critical points of the gauge 
model) and should abruptly break down for {3 < (3 r . It is important to recall that in the 
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Table 6: Results of the fit according to eq. In the first three columns the values of 

(3 and L of the data used in the fit and the corresponding value of T/T c are provided. In 
the fourth column, we give the range of values of R included in the fit. In the fifth column 
the reduced x 2 is given. The last two columns contain the best fit results for a and for the 
Luscher term. 



R-min 




o 


c 


8 


181 


0.0105622(7) 


0.824(1) 


10 


47 


0.0105454(8) 


0.883(2) 


12 


15 


0.0105364(9) 


0.922(2) 


14 


6.6 


0.0105315(10) 


0.949(3) 


16 


3.9 


0.0105285(11) 


0.967(5) 


18 


1.2 


0.0105255(12) 


0.990(6) 


20 


1.1 


0.0105247(14) 


0.996(8) 


22 


1.0 


0.0105241(15) 


1.003(10) 



Table 7: Results of the fit according to eq. lfI3)) for the data at (3 = 0.75180, L = 80. The 
first column contains the minimum value of R included in the fit (recall that for this value 
of (3 the correlation length is twice the one of f3 = 0.73107 and a similar doubling of all 
the length scales must be kept into account when comparing the data at the two values of 
(3). In the second column the reduced x 2 is given. In the last two columns we provide the 
best fit results for a and for the coefficient in front of the bosonic string correction. 

above fits we neglected the short distance data. As we shall see below the short distance 
deviations from the free string behaviour are instead affected by scaling corrections. 

In the following section we shall make a further step in our analysis. We shall perform 
some kind of "absolute test" of the effective string, combining the data so as to eliminate 
the string tension from the game as well, thus only leaving the effective string corrections 
in the data. 

2.2.3 An "absolute" test of the effective string picture 

This can be easily done by building the combinations 

H(R U R 2 ) = -- log [ G{Ri)G{R2 + 1} ) (14) 
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b=0.75180, L=80 
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Figure 1: Comparison between our values of H(R, 20), the bosonic string prediction, eq. (0) 
(continuous line) and the 1/R approximation to eq. © (dashed line) for the data at 
P = 0.75180, L = 80. 



in which the string tension disappears. We plot in fig. ^ and El a sample of the results we 
obtain. In particular we concentrate on the sample at the largest value of (3 and fix one of 
the two entries to the first acceptable value of R m in as extracted from the previous analysis, 
i.e. i?2 — 20. In fig.s ^ an d 12] we compare our results with the bosonic string prediction 
for H(R, 20) obtained using eq. © and also with what one would obtain (dashed lines in 
fig.s^andEJ) keeping only the 1/R correction in eq. ©. We find for R\ > 20 an impressive 
agreement between our data and the effective string predictions. All the data agree within 
the errors. Instead, for R\ < 20 the large deviations from the free effective string appear. 
We shall discuss these deviations in the forthcoming sections. It is very interesting to 
observe that, in agreement with the results of the first level analysis discussed above, 
for large enough values of R the contribution due to the next to leading, exponentially 
decreasing, terms of eq. © is clearly visible in the data. It is clear that these terms cannot 
be neglected in the fits. 
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b=0.75180, L=80 




Finally in order to have a global view of all our data we built the differences 

A(R h R 2 ) = H(R h Ra) - i (F^R, + 1,L) + F^(R 2 , L) - F^(R 2 + 1, L) - F^(R U L)) . 

(15) 

With this definition, we have A = for all choices of R\ and R 2 if the free effective string 
picture of eq. (jOJ) is the correct description of the data. We report in fig. El as an example 
of the results that we obtained the histogram of A for (3 = 0.75180, L = 80 for all the 
values of Ri,R 2 > 20. The overall agreement between the data and the effective string 
prediction is very good. With the only exception of the L = 20 set the data show a 
rather nice Gaussian like distribution around the expected A = value. The variance of 
these Gaussians is of the same order of the errors of H(Ri,R 2 ) extrapolated from those 
reported in tab. El and El and no systematic deviation or asymmetry is visible in the data. 
These observations suggest that the (very small) spread around the A = value is a purely 
statistical effect (as a side remark this is also a test of the fact that our statistical errors are 
reliable). On the contrary the sample at L = 20 shows some systematic asymmetry. This 
is an indication that for T > T c /5 the large distance, high temperature effects discussed 
in jlj start to give non-negligible contributions within the precision of our data. 

2.2.4 The zero temperature string tension 

As a byproduct of the previous analysis we may obtain from our fits a very precise deter- 
mination of the zero temperature string tension for three sets of data. They are reported 
in tab. [HI where we have listed the results of the fits according to eq. (fT5Jl . 



p 


L 


T/T c 


Rmin 




a 


c 


0.65608 


20 


1/10 


9 


0.58 


0.204864(9) 


1.017(11) 


0.73107 


80 


1/20 


10 


1.3 


0.044023(3) 


1.003(4) 


0.73107 


40 


1/10 


10 


2.0 


0.044019(4) 


1.005(5) 


0.73107 


20 


1/5 


10 


0.67 


0.043985(4) 


0.988(6) 


0.75180 


80 


1/10 


22 


1.0 


0.0105241(15) 


1.003(10) 



Table 8: Results of the fit according to eq. J23J). In the first three column the values of 
(3 and L of the data used in the fit and the corresponding value of T/T c are summarized. 
In the fourth column the minimum value of R included in the ht is quoted. In the fifth 
column we give the reduced y 2 . The last two columns contain the best ht results for a and 
for the coefficient in front of the effective string correction. 

The values of a that we find for (3 = 0.75180 is slightly smaller than the value we 
used in our previous papers [3111, which was a = 0.010560(18). The effect is very small, 
but nevertheless it is outside the quoted error bars (while in the case of (3 — 0.73107 the 
error bars are too large to detect this effect). An interesting consequence of this (very 
small) rescaling of a is that it could explain the residual (small) systematic deviation we 
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Figure 3: Histogram of A(Ri, R 2 ) for all the values of R±, R 2 > 20. Notice the very small 
scale of the A axis. 
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found in {4J between our data and the Nambu-Goto effective string truncated at the first 
perturbative order (see fig.s 2, 3, 4 and 5 of jl]). Both the sign of these deviations and the 
fact that they are proportional to L are in agreement with this explanation. 



3 Higher order corrections at short distance 

With the above results at hand, we are in the position to address the issue of the higher 
order correction at short distance in a reliable and precise way. We shall perform this anal- 
ysis in two steps. First we shall study the samples at (3 = 0.73107 and (3 = 0.75180 trying 
to quantify the higher order terms in the potential, without referring to any particular 
effective string picture. Then, in the second step we shall compare our results with the 
expected behaviour of the Nambu-Goto string. In this second step we shall also include 
the data at j3 = 0.65608 in our analysis. 



3.1 Higher order corrections 

We decided to study the higher order corrections by looking at the observable (see the def. 
of eq. (HHJ)) A(R,R-n). 

This quantity is interesting for several reasons. First of all, as we mentioned above, it 
is different from zero only if higher order corrections are present. Moreover, for n = 1 we 
may easily relate it with the observable c(f) introduced in j2J. The relation is 

H(R,R-1) = 2^- (16) 

(for a definition of f see As a side remark, this means that H(R,R — 1) (suitably 
normalized) is an estimator of the central charge of the underlying conformal field theory 
(CFT), (see [4 for a discussion of the CFT description of effective string theory). Thus 
A(R, R — l) measures the deviation of the central charge of this CFT from the free bosonic 
value (i.e. c = 1). 

In order to compare values of this quantity for samples at different values of a and for 
different choices of n it is useful to define the following scale invariant version of A 

In building this observable we used the high precision results for a obtained in the previous 
section. 

We are primarily interested in identifying the exponent of the correction, thus we fitted 
D(R) with the ansatz 

D i R ) = fS§ ■ ( 18 ) 
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It is possible to relate g a +2 to a corresponding higher order correction in the free energy. 
More precisely, if the free energy has a correction of the type 



then in D(R) we expect to find a correction proportional to g a+ 2/R a+2 (plus higher order 
terms due to the lattice discretization) with g a+2 = oc(a + l)j a - Notice that in eq. (fT9*)l 
we have rescaled the coefficient of the higher order correction with a suitable power of a 
so as to make it adimensional. We also implicitly assumed that the higher order term 
is proportional to L. We shall see below that this assumption agrees with the numerical 
results. 

Looking at the discussion of sect. ^ it is easy then to relate the presence of a nonzero 
7 Q term to a precise contribution in the string action. In particular, if the fit with a = 1 
gave an acceptable (i.e. order unity) reduced x 2 and a value 71 7^ 0, then this would 
indicate that the Liischer term does not properly describe the data and that a different 
coefficient in front of the 1/R correction to the interquark potential is needed in this short 
distance regime. If instead the a = 2 fit gave an acceptable Xri then this would indicate 
that a boundary term might be present in the data, and we could estimate the boundary 
parameter b (we shall discuss the boundary contribution in sect. 0] below). Notice that, 
according to our normalization choice, the parameter b introduced by Liischer and Weisz 
in [2] is given by b = — ^ - . We shall discuss this issue in detail below. Finally, 

an acceptable xl with a = 3 would indicate that higher order terms are present in the 
effective string action. In particular, the Nambu-Goto action would precisely give such 
a contribution (see below). The value of 73 extracted from the fit could help to identify 
the form of these higher order terms. Clearly, once for a given a we find 7 a 7^ all the 
higher powers of R in the fit would have nonzero coefficients, regardless of the presence of a 
corresponding term in the action, simply because of the lattice artifacts. It is well possible 
that a mixture of all these terms is present in the action, but unfortunately our data are 
not precise enough to allow for more than one parameter fits. Thus we shall only be able 
to identify (if it exists) the leading term among them. As a consequence one should look at 
the results which we shall now discuss more as a qualitative indication than a quantitative 
estimate of the higher order terms in the string action. Notwithstanding this, a few precise 
pieces of information can be obtained from our data. 

As a last point let us mention that it is important to properly select the range of values 
of R to be included in the fit. According to [3] and jlj, we must choose R > R c (see the 
values reported in tab. HJ). 

We performed the fit for various integer values of a; the results are listed in tab. IDlfTTI 
(see also fig. |U and fig. EJ) for all the data except those at L = 20 (which, as we have shown 
in the previous section, seem to be affected by high temperature corrections). 

Let us briefly comment on these results. 

1] The three sets of data show a remarkably similar pattern. 




(19) 
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Table 9: Results of the fits according to eq. (JTSJ) for the data at (3 = 0.75180 and L = 80. 
In the first column we give the value of a, in the second column the reduced x 2 , in the last 
two columns the values of g a+2 an d j a - 



a 




9a+2 


la 


1 


2.8 


0.033(3) 


0.0166(13) 


2 


1.8 


0.053(4) 


0.0089(7) 


3 


1.5 


0.079(6) 


0.0066(5) 


4 


1.5 


0.114(9) 


0.0057(4) 



Table 10: Same as tab. Elbut for the data at (3 — 0.73107 and L = 80. 



a 


"X,r 


9a+2 


la 


1 


22.1 


0.037(1) 


0.0182(4) 


2 


6.8 


0.052(1) 


0.0087(2) 


3 


1.5 


0.070(1) 


0.0059(1) 


4 


1.0 


0.092(2) 


0.0046(1) 



Table 11: Same as tab. El but for the data at (3 = 0.73107 and L = 40. 



a 




9a+2 


la 


1 


10.4 


0.036(1) 


0.0180(6) 


2 


4.1 


0.052(2) 


0.0087(3) 


3 


1.8 


0.070(2) 


0.0058(2) 


4 


1.4 


0.092(3) 


0.0046(1) 
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2] Looking at the \ 2 values it is clear that our data exclude a a = 1 correction. This means 
that the Liischer term we have subtracted perfectly describes the 1/R behaviour of 
the data. 

3] A boundary term (a = 2) gives reduced x 2 's which are higher than the a > 2 ones in all 
three cases. If we assume that the data are described by a boundary correction, then 
the amplitude of such a correction turns out to be of the same order of magnitude 
(and sign) of those that we find in the large distance regime (see the discussion at 
the end of sect. EJ), but the scaling behaviour of the large distance results is not the 
one expected for a boundary term. 

4] The fits seem to support a value between 3 and 4 for a. This result is compatible 
with a 1/R 3 "Nambu-Goto like" correction or, more likely, a mixture of a 1/R 3 plus 
higher order corrections. 

5] The two sets of results at (3 = 0.73107 are compatible within the errors. This means 
that also the higher order corrections in the free energy F(R, L) scale linearly with 
L. 

6] Comparing the sets of data at (3 = 0.73107 and the one at (3 = 0.75180 we see a good 
scaling behaviour in (3. This is also evident if one looks at fig. 0] in which the two 
sets of data are plotted together. 

3.2 Comparison with the Nambu— Goto string 

The most impressive feature of the above analysis is that, although the power that we find 
1/R? is exactly the one expected according to the Nambu-Goto action, the value of the 
coefficient 73 ~ 0.006 that we find is definitely different with respect to the Nambu-Goto 
prediction eq. ©, namely 73 = = ~0. 00856 , which is opposite in sign and more or 

less 4/3 in magnitude with respect to our fit's results 3 . 

In fig.sEJandElwe plot our data for the three values of (3 (for (3 = 0.73107 we chose the 
sample with L = 40). For comparison, we also plot the Nambu-Goto expected behaviour 
(solid line) and the two curves obtained using the two best fit values 7 = 0.0059 and 
7 = 0.0066 for (3 = 0.73107 and (3 = 0.75180 respectively (dashed lines). 

Looking at these figures, we see that the data for /3 = 0.65608 show a large scaling 
violation. This had to be expected due to the rather small value of the correlation length 
£ < 1 for this (3. Note that even the free energy of a free field theory on a lattice shows 
corrections of the order 1/R 3 . As a little exercise, we have computed the free energy of the 
free field theory on a square lattice with the appropriate boundary conditions. We have 
used the most naive discretisation of the derivate of the field. We would expect that the 
artifacts in the Ising model are at least as large as those in the free field theory on the 

3 A similar disagreement in the short distance regime was recently reported in [§]. 
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lattice. Indeed we find that the corrections in the free field theory have the same sign as 
those of the Ising gauge model at (3 = 0.65608, but a smaller amplitude. 

On the other hand, the fact that the deviations between (3 = 0.73107 and (3 = 0.75180 
are much smaller than those with (3 = 0.65608, makes us confident that for (3 > 0.73 
lattice artifacts have only a minor impact on the higher order corrections and hence we 
can extract continuum physics results. 

A possible interpretation of this scenario, keeping also into account the results of our 
previous paper [3] is that, while the Nambu-Goto effective string provides a correct de- 
scription of finite temperature effects of the interquark potential at large distances, while 
it fails to describe corrections at small distance R. Along with Dirichlet boundary condi- 
tions, some irrelevant operator becomes important and drastically modifies the shape of the 
interquark potential (or better: the shape of its deviation from the free string behaviour). 

The simplest candidate for such a non-universal behaviour is a boundary term of the 
type discussed in the next section. A positive and large enough value of b could perhaps 
justify the observed behaviour of D(R). However as we shall see in the next section the 
large distance behaviour of the interquark potential in the Ising case seems to suggest 
that b ~ 0. Thus it is possible that some other operator, maybe an extrinsic curvature 
contribution [20, 21J is actually responsible for the observed behaviour. In order to better 
understand this point we shall devote the following two sections to a detailed discussion of 
the boundary term. 

An alternative scenario is that the short distance breakdown of the Nambu-Goto string 
description has to be explained by effects beyond an effective string theory; for a discussion 
we refer the reader to ref. 0. 

4 The boundary term 

As we have seen in the previous section, in order to understand the short distance be- 
haviour of the effective string potential it would be important to have an independent 
estimate of the boundary parameter b. This can be achieved by studying the large dis- 
tance data. However, in order to perform this analysis, we must know the functional form 
of the boundary correction in the large distance regime. This section is devoted to such 
calculation, while in the next section we shall compare this result with the large distance 
data of ref. 0j. 

4.1 The functional form of the boundary correction 

As mentioned in the introduction, the most general effective string action in presence of 
boundaries requires the inclusion in the action of terms that are localized at the boundary. 
The simplest possible term of this type is [2j 
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Comparison of three samples 
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Figure 4: Values of -D(-R) versus R^/a for the sample at f3 — 0.65608 (white squares), 
the one at (3 — 0.75180 (black squares) and the L = 40 one at j3 — 0.73107 (crosses). 
The continuous line is the Nambu-Goto expectation, the two dashed lines are the best fit 
results for (3 = 0.73107 and (3 = 0.75180 discussed in the text. 
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Comparison of three samples 
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Figure 5: Same as fig. HJbut with a higher resolution. 
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Ab = \ I dt 
4 Jo 



dhy fdh 

dz) z=0 \dz y ,, 



(20) 



where b is a parameter with dimensions [length], h denotes the transverse displacement of 
the string (see the appendix for a detailed discussion) and the \ factor has been added to 
agree with the conventions of 

This additional term can be treated in the framework of the zeta-function regularization 
in a way similar to the Nambu-Goto one. We report in the appendix a detailed account 
of this calculation and only list here the final result, which turns out to be rather simple. 

At first perturbative order in b, the regularization of the free string action plus the 
boundary term gives the same result as the pure free string action, namely the Dedekind 
function discussed above, provided one replaces the interquark distance R with 

R^R* = — -7— y (21) 

(remember that b has the dimension of a length). 

Thus, if we denote the effective string contribution to the interquark potential at this 
order (with an obvious choice of notation) as Fg(R, L) + F b (R, L), we find (we shall keep 
from now on d = 3 to simplify the equations) 



F}(R,L)+F b q (R,L) = logrj(i^j . (22) 



Expanding this result in the short distance regime (i.e. for 2R* < L), and neglecting the 
exponentially decreasing corrections we find 

F'(R,L) + F b q (R,L) = -^(l + 2^y , (23) 
which becomes at the first order in b 

F q 1 (R,L) + F b q (R,L) = -^(l + ^) , (24) 
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which is exactly Liischer and Weisz's result 0. 



4.2 Large distance behaviour of the boundary correction 

The major advantage of having the complete functional form of the boundary correction is 
that we can now look at its large distance behaviour. By using the modular transformation 
of the Dedekind function we find for 2R* > L 
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F'(R,L) + F b (R,L) 



ttR* 



1 9 7?* 

~" ^ n=l 



(25) 



Neglecting again the exponentially decreasing terms, and keeping only the first order in b 
we end up with 



f}(r,l)+f;(r,l) 



TlR 
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H) 



1 

2 
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IT 
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1 + 2- 
R 



(26) 



which means (remember again that we only keep the terms proportional to b in the 
expansion) 



(27) 



from which we obtain 



nb 



F(R,L) = +—- 



b 

6L 2R 

Thus we see that the effect of the boundary term is to renormalize the constant term 
k(L) (which however is not relevant in the ratios evaluated in Monte Carlo simulations) 
and to add a further 1/R correction. 

A few observations are in order at this point. 

• This result is rather satisfactory from a physical point of view: the effect of a bound- 
ary term can only decrease as a function of R (i.e. when the two boundaries are 
far apart), and this is indeed the case in the large R regime, too. This is in sharp 
contrast with the behaviour of the string fluctuations, whose dominant term in the 
large distance regime is always (i.e. at any order in the 1/aRL expansion) a linearly 
rising correction. 

• Although both in the short distance and in the large distance regimes the correction 
induced by the boundary term has the same sign, nevertheless there is a well defined 
change of behaviour: from — ^f^ to — This feature may simplify the identification 
of such a term in numerical simulations. Notice in particular the lack of L dependence 
in the large R regime. 



5 Comparison with the Ising gauge model in the large 
distance regime 

The aim of this section is to compare the prediction of eq. ()27j) with the results of the 
simulations reported in [3] for the large distance regime of the interquark potential in the 3d 
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gauge Ising model. In particular, we shall concentrate on the set of simulations performed 
at (3 = 0.75180, which were the most precise ones. The simulations were performed on 
lattices with space-like size N s = 128 and "time" sizes L = 10, 12, 16, 24, corresponding to 
|T C , |T c , \T c and |T C , respectively. For each temperature, we studied the correlators for 
various values of the interquark distance R, ranging from R = 8 to R = 48. The data are 
reported in 

The most effective observable to identify the possible presence of a boundary correction 
is the following combination 



W) = - log g( * * 1} ~ Fi(R + 1) + Fi{R) 



. G(R + 1) , 
l ° g -G(RT + hg 



L \ 



lot 



V 



2(R+1)J 



(29) 



In order to understand the results of the fit it may be useful to extract from eq.s (@J) 
and (jHJ) the expected behaviour in the large distance regime for this quantity (which as 
we mentioned several times is drastically different from the short distance one discussed in 
sect.sElandEJ). Discarding exponentially decreasing terms, we find 

F}{R + 1, L) - Fg(R, L) = + \ log (30) 

and (see eq.s (2.37) of ref . [I]) 



TV 2 



F„(R + 1, L) — F„(R, L) — ; -. (31) 

qK ' ; qK ' ' 72aL 3 8aLR(R+l) V ; 

These corrections should be compared with what is expected from the boundary term, 
namely 

Fq(R + 1,L)- F b q (R, L) = 2R{ x +1) ■ (32) 

The boundary correction to Q2{R) has the same R dependence as the Nambu-Goto term, 
but a different L dependence (this is an obvious consequence of the fact that b is a di- 
mensional parameter): that makes it very easy to disentangle between the two corrections. 
In particular, fitting them with the data of 0], we find that b is compatible with zero, i.e. 
that no boundary correction seems to be present in the 3d Ising gauge model. 

More precisely, we fitted Q2(R), for each of the four values of L, with the following 
formula 

g 2 ( J R) =ao + &o ^ 1 +1) , (33) 

where ao is related to the string tension, while bo encodes the large distance 1/R correction 
to the potential we are interested in. With this choice of normalization, if the 1/R correction 
is completely due to the boundary term, then b = 26 - Since for this value of j3 we have 



25 



Table 12: Results of the fits to eq. In the first column the value of L, in the 

second column the reduced \ 2 JS given, in the last two columns the values of bo and 
the corresponding expectation — see eq. — according to the Nambu-Goto string. 



L 




bo 


1/BaL 


10 


2.72 


-1.31(7) 


-1.184 


12 


2.33 


-0.97(6) 


-0.986 


16 


2.30 


-0.62(5) 


-0.740 


24 


2.54 


-0.30(4) 


-0.493 



R c ~ 12, we only fitted the data for R > 16, for a total of 6 points. The results of the fits 
are reported in tab. \H\ 

Few observation can be made on these results 

a] The magnitude of the reduced v 2 (which is remarkably constant as L changes) indicates 

that in the large distance regime the Nambu-Goto string (truncated at the first 
perturbative order) is a reasonable approximation, but does not fully describe the 
data. 

b] There is a clear L dependence in the best fit values of &o- These values qualitatively 

agree with the Nambu-Goto expectation, but they seem to show an even steeper 
dependence on L. Apparently little room is left for a boundary correction term. 
A naive extrapolation to L — > oo suggests a very small value of b, most probably 
compatible with zero. 

c] It is interesting to compare what we find here with the fits discussed in sect. El (see 

tab.sEl EH and ED we find for this value of (3: b = b/2 = -f ^ ~ -0.35 which 
is of the same order of magnitude and has the same sign of the values reported in 
tab.EEl 

6 Concluding remarks 

Let us summarize the main results of our analysis. 

a] We obtained an explicit form for the boundary correction, at first perturbative order 

in b, which is valid for any value of R and L. The importance of this result is that 
in the large R regime, by virtue of the peculiar form of the boundary correction, i.e. 
the fact that it is L-independent, it is much simpler to measure the value of b. 

b] As a first application of this result, we evaluated b in the 3d gauge Ising model, using the 

large distance data published in It turns out that b is very small and compatible 
with zero. 
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c] In order to test the effective string picture in the short distance regime as well, we 

performed a new set of simulations for the 3d gauge Ising model. A careful analysis 
of the data in the range of T and R values in which higher order terms can be 
neglected shows a perfect agreement with the prediction of the free bosonic string. 
However it turns out that in order to describe the data one must keep into account 
the whole functional form of the correction. Cutting it to the Liischer term 
only may induce a relevant mismatch with the data. 

d] Finally we studied the higher order corrections for small values of R. In this regime 

the Nambu-Goto correction, truncated at the first perturbative order, should behave 
as 1/R 3 . Assuming 6 = and fitting the data with this correction (the a = 3 fit in 
the sect. EJ) we found an acceptable Xri but the wrong coefficient. 

This is indeed rather a puzzling result since in the large distance, high temperature 
regime we found instead quite a good agreement with the Nambu-Goto prediction (see 
[3]). In principle, a description holding at short distances could even be a picture different 
from an effective string However, it is also possible that the short distance breakdown 
of the Nambu-Goto string scenario could be explained within the framework of an effective 
string theory. In particular, following the comments that we made at the end of sect. El 
we see two possible explanations for this behaviour. 

1] The Nambu-Goto picture could be the wrong assumption. In this framework the 
good agreement at large distance is only a coincidence, and one should look for 
some other, more exotic, effective string action to simultaneously describe the short 
distance as well as the finite temperature behaviour of the potential. Notice however 
that none of the generalized strings discussed in ^2] satisfies these constraints, thus 
some drastically different model would be needed to follow this line. 

2] The short distance deviation with respect to the Nambu-Goto prediction might be 
due to an irrelevant operator which may possibly have a string-like description. Two 
natural options are: 

• a higher order boundary term which decreases with a higher power of R and 
could thus be compatible both with the behaviour that we extract from the 
short distance fits (which suggests a power higher than 1/R 2 ) and with the 
large distance regime. 

• a term proportional to the extrinsic curvature of the string |2*U1 |2*T] . It is pos- 
sible that such a term gives significantly different contributions in presence of 
Dirichlet and periodic boundary conditions. 

In this framework it could be useful to compare the short distance behaviour of 
different models to see if some common feature emerges. The nice feature of this last 
option is that it does not require a large value of the b parameter and thus it does 
agree with all our numerical results. 
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Following the suggestion of the last point above we tried to test if some common 
behaviour is shared by different gauge theories. To this end, we performed the same 
analysis discussed in sect. El with the data on the 3d SU(3) gauge model published in j2]. 
In particular, we concentrated on the sample at ft — 20, see the data in tab. 3 of ref. |2], 
which correspond to a value a ~ 0.0342. In agreement with the results of |2|, we found 
that the data are well described by the free string correction only, with apparently no 
need of higher order corrections. The most impressive signature of this behaviour is given 
by the observable D(R) introduced in sect. HU which is different from zero only if higher 
order corrections are present. It turns out that for SU(3), D(R) is almost compatible 
with zero within the errors, for all of the values -JaR > 1.3. To complete the comparison 
we also studied a set of data obtained in the 3d SU(2) model taken from [221 • These 
data are characterized by a slightly smaller string tension a ~ 0.0263 and a temperature 
T/T c ~ 1/7 (see [22| for details). The value of D(R) for these two sets of data, together 
with the L = 80, ft = 0.73107 sample of the Ising model are reported in fig.s and 
Notice that the SU(3) and SU(2) samples have values of a which are smaller than the Ising 
one. Looking at these figures one sees that the SU(2) data lie somehow in between the 
Ising and the SU(3) ones. This is also confirmed by the 73 value extracted from the SU(2) 
data: 73 ~ 0.004 [221 which is smaller (and outside the error bars) than the one that we 
obtained in the present paper for the Ising case. 

A similar comparison between the 3D SU(2) and Z 2 gauge models for the first excited 
string level can be found in ref. 0. In this case, the SU(2) and Z 2 data seem to be 
compatible. 

Further studies are needed to better characterize these behaviours, to understand if 
they may be related to some feature of the gauge group (say, for instance, to the center of 
the group as it is suggested by the analogy between Ising and SU(2) behaviours) and to 
see if they support any of the scenarios proposed above. 

Acknowledgements. This work was partially supported by the European Commission 
TMR programme HPRN-CT-2002-00325 (EUCLID). M. P. acknowledges support received 
from Enterprise Ireland under the Basic Research Programme. M. H. thanks PPARC for 
support under grant PPA/G/O/2002/00468. 



28 



Comparison of SU(2), SU(3) and Ising 
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Figure 6: Values of D(R) versus R\fo for the d=3 SU(3) gauge theory (white squares), the 
Ising gauge model (black squares) and the 3d SU(2) gauge theory (crosses). The vertical 
line corresponds to R — R c . 
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Comparison of SU(2), SU(3) and Ising 
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Figure 7: Same as fig. El but with a higher resolution. 
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A Evaluation of the functional form of the boundary 
correction. 



Let us begin by considering the case when we assume that the dynamics of the effective 
string world sheet is described by a pure Nambu-Goto action, neglecting a possible bound- 
ary term: under this hypothesis, the action is given by the string world sheet surface area 
(measured in units of the inverse string tension), and the expectation value for the Polyakov 
loop correlator can be expressed as the partition function: 



Z 



[Vh}e 



-aA 



(a.i; 



The Nambu-Goto action Aq can be written in the following form: 



pR pL 

Aq= dz dt y/det g a p 
Jo Jo 



(A.2) 



where g a n is defined as: 



g af3 = d a X*dpXp (A.3) 
Let us focus on the case of three dimensions: we can parametrize the X field in the following 



way: 



X x 

x 2 

X 3 



t; 
h 



h(z,t) ; 



(A.4) 



(A.5) 



and the h field is associated to the string world sheet surface transverse displacement with 
respect to a minimal area plane. Thus we have: 



An 



dz f dt\ll + (^] + (^r 
o Jo 



\dz J 



dt 



(A.6) 



A Taylor expansion of the square root appearing in Eq. (jA.6|) gives — to first non-trivial 
order: 



[jh) + \ dt) 



and the partition function can be approximated as: 



Zo+i 



-aRL 



[Dh] e 



./ ^/>[(f) 2 + (t) s 



(A.7) 



(A. 
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On the other hand, if we allow for the existence of a possible "boundary term" , the effective 
action will take the form: 

r [Vh] e- aA (A.9) 



with: 



A = A + - 



L 

dt 



o 



fdh\ 



fdh\ 



{dz) z=0 + {dz) 



z=R 



(A.10) 



where b is a parameter with dimensions [length], while the \ factor is just a conventional 
choice (see also j2]). 

Let us notice that, in the case of the system we are presently interested in, namely: a finite 
temperature confining gauge theory 4 , the ^-dependence of the eigenfunctions appropriate 
to the case of the pure Nambu-Goto action can be factored out as: 



thus: 
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dh 

dz 



( 



nnz\ 



V R 



-R CO \ R 



, n G N 
/mrz\ 



(All) 
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(A.12) 



where we omitted the eigenfunction term expressing the t-dependence; thus we have: 

2 / or \ 2 
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fdh\ 

n 2 ix 2 



n 2 7c 2 



z=R 



R 2 



cos 



'nirz\ 



V R 



1 



2 = 



n 2 7T 2 



cos 



'nirz\ 



V R 



z=R 



(A.13) 



As it is concerned with the momentum in the compactified direction, the eigenvalues are: 

Air 2 m 2 



L 2 



If we define: 



P 

T 



m e Z 



R 

L 
2 



the eigenvalues of the (—A) operator can be written as: 

,2' 



7l~ 



n m 

-7 + 



(A.14) 
(A.15) 

(A.16) 



.P T 

4 Thc string world sheet associated to the expectation value of the Polyakov loop correlator in a finite 
temperature gauge theory is characterized by fixed boundary conditions in the z direction, and periodic 
boundary conditions in the compactified, t direction. 
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Taking into account the boundary term too, the complete action A reads, to first non-trivial 
order: 



A~ dz l dtll , 
Jo Jo I 2 



dhy fdh\ 
dz) + \dt) 



b f 
+- / dt 



dh 

dz 



+ 

2 = 



dh 

dz 



(A.17) 



Now, we propose to work out the correction to the V(R) interquark potential to first order 
in b in a perturbative expansion; to do this, we have to evaluate: 



J o+i 



[Dh] e- a - 



A 



(A.18) 



J up to first order 



The Taylor expansion of the exponent gives: 

a f R f L 
a A ~ oRL + - dz dt 



o 

ab f L 7 
+ — / dt 

4 Jo 



fdh\ 2 fdhX 
\dz~) + \di) 



(dhY 



dhV 

dz )z=R 



(A.19) 



We can think about this formula as the starting point in a perturbative analysis in terms 
of the b parameter, namely we can consider the boundary term as a perturbation of the 
pure Nambu-Goto action Eq. (|A.2|) . 
First of all, let us notice that: 

1. by virtue of the fact that the eigenvalues of the unperturbed problem are known - 
see Eq. (IA.16|) — the double integral appearing in Eq. (|A.19|) can be written in a 
simpler way: 



7T 



R [-L 

dz I dt 
o Jo 



1 1 

n / n m 
1 

p2 r 2 



fdh\ 2 (dh 



R pL 

dz / dt{...) [zt? 



(A.20) 



where the notation (. . Ofz.tp represents (the square modulus of) the complete, unper- 
turbed eigenfunction. 



2. On the other hand, Eq. (jA.13|) shows that the two addends appearing in the integral 
associated to the boundary term in Eq. (|A.19|) are equal, and they give rise to a 
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contribution which can be written as: 
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dt 



Oh 
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dz 1 n \dz 
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dt 



(A.21) 



In this case, the integrand 
part of the eigenfunction 5 . 



is just (the square modulus of) the t-dependent 



Thus, in the "pure Nambu-Goto" case, aAo involves terms like: 
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dz 



L 

dt (. . .) 
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but when the boundary term is included this expression has to be replaced by: 
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ah n 2 TT 2 
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The integrals appearing in Eq. ()A.23jl give rise to results 6 in the form: 

C 2 -RL 
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dz 







L 

dt(..:, 
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dt\ 



o 



[z,t] 2 



where C 2 and G\ are just pure numbers 7 ; thus we can rewrite Eq. ()A.23jl as: 
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Defining: 
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Ci 
C 2 



(A.22) 



(A.23) 



(A.24) 
(A.25) 



(A.26) 



(A.27) 



5 The z-dependence does not appear, because we already evaluated its value in z = and in z = R. 

6 The results appearing in Eq. ljA.24f> and Eq. (|A.25|) refer to an unnormalized eigenfunction, but this 
is not relevant to our present calculation, since the only difference in the normalized eigenfunction case is 
that both quantities are multiplied by a common factor, and their ratio is unchanged. 

7 Strictly speaking, this is true only if we do not consider the usual normalization factor; however, even 
in the more general case, the ratio between Ci and C\ will be the same as we find here. 
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then Eq. (lA~26j) can be rewritten as: 

7T 2 

a— ■ C 2 RL 
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7T 
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(A.28) 



By defining: 



Eq. ()A.28|) takes the form: 
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(A.29) 



(A.30) 



In the pure Nambu-Goto theory (with no boundary term), the corresponding term: 



a— ■ C 2 RL 
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2 2 

m n 

^ 2 ~ + R* 



(A.31) 



is associated to the Z\ term in the partition function factorization formula Eq. (|A.8|) : in 
three dimensions, Z\ turns out to be equal to: 



Z x 
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^ l ^R 



(A.32) 



and its (additive) contribution to the interquark potential is: 

V 1 (R) = ~\ogZ 1 = ^]ogri(i^ 



L 

Since Dedekind's function is defined as: 
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(A.33) 



(A.34) 



n=l 



In the regime ^ ^> 1, we can approximate: 
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(A.35) 
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which is nothing but the celebrated "Liischer term" 
When the boundary term is included in the action, we get a similar result, provided we 
replace: 

R — ► R* = R (l + k^j 2 (A.36) 

We can evaluate k in simple way: its definition is given by Eq. (|A.27|) . which involves the 
integrals in Eq. (jA.24j) and in Eq. (jA.25j) . thus: 

A; '! (a - 37 » 

Since the two eigenfunction parts expressing the t-dependence and ^-dependence are fac- 
torized, the denominator integrals factorize, too. Moreover, we know that the z-dependent 
part of the eigenfunction is expressed by Eq. (|A.lljl ; these observations allow one to write: 

k = R R J ° dt{ ---y = R. 1 L = 2 (A.38) 

where we exploited the fact that n G No- 

Thus, in the case when we allow for the presence of the boundary term, the first-order 
non-trivial contribution to the interquark potential will be: 

V^hogrjU^j (A.39) 

It is important to stress the fact that this result does not rely on the assumption to be 
in the "large ^ regime": as a matter of fact, the only hypothesis for the validity of Eq. 
flA.36|) is that h can be considered as a perturbatively small parameter. 
On the other hand, if we restrict to the case when ^ ^> 1, then Eq. (jA.39|) can be 
approximated by: 

V* : ^ 



2AR+ 



24i?(l + 2A)- 



" 1 + 4^-^(1 + 4) (A.40) 



2AR V R 2AR \ R 
which exactly reproduces the term appearing in Eq. (3.9) in j2]. 
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